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Abstract. We use simultaneous observations of tropospheric 
ozone and outgoing longwave radiation (OLR) sensitivity to 
tropospheric ozone from the Tropospheric Emission Spec- 
trometer (TES) to evaluate model tropospheric ozone and 
its effect on OLR simulated by a suite of chemistry-climate 
models that participated in the Atmospheric Chemistry and 
Climate Model Intercomparison Project (ACCMIP). The en- 
semble mean of ACCMIP models show a persistent but mod- 
est tropospheric ozone low bias (5-20 ppb) in the South- 


ern Hemisphere (SH) and modest high bias (5-10 ppb) in 
the Northern Hemisphere (NH) relative to TES ozone for 
2005-2010. These ozone biases have a significant impact on 
the OLR. Using TES instantaneous radiative kernels (IRK), 
we show that the ACCMIP ensemble mean tropospheric 
ozone low bias leads up to 120 mW m -2 OLR high bias lo- 
cally but zonally compensating errors reduce the global OLR 
high bias to 39 d= 41 m Wm -2 relative to TES data. We show 
that there is a correlation ( R 2 = 0.59) between the magnitude 
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of the ACCMIP OLR bias and the deviation of the ACCMIP 
preindustrial to present day (1750-2010) ozone radiative 
forcing (RF) from the ensemble ozone RF mean. However, 
this correlation is driven primarily by models whose absolute 
OLR bias from tropospheric ozone exceeds lOOmWm -2 . 
Removing these models leads to a mean ozone radiative forc- 
ing of 394 =b 42 m Wm -2 . The mean is about the same and 
the standard deviation is about 30 % lower than an ensemble 
ozone RF of 384ib60mWm -2 derived from 14 of the 16 
ACCMIP models reported in a companion ACCMIP study. 
These results point towards a profitable direction of combin- 
ing satellite observations and chemistry-climate model sim- 
ulations to reduce uncertainty in ozone radiative forcing. 


1 Introduction 

Tropospheric ozone plays a central role in both atmospheric 
chemistry and climate. Due to significant increases in an- 
thropogenic emissions of its precursors since preindustrial 
times, tropospheric ozone has the third highest impact as an 
anthropogenic greenhouse gas in terms of radiative forcing, 
350 [250-650] mWm~ 1 2 (Forster et al., 2007), but is distin- 
guished from other greenhouse gases by having significant 
spatial and temporal heterogeneity due to its relatively short 
lifetime. Ozone is an essential part of the oxidative capac- 
ity of the troposphere, both as an oxidant in its own right, 
but especially as a precursor of the hydroxyl radical (OH) 
(Wang and Jacob, 1998; Voulgarakis et al., 2013). Reaction 
with OH controls the chemical lifetime of many species in 
the atmosphere, including methane, providing a link between 
ozone concentrations and the global warming potential of 
methane (Shindell et al., 2005, 2009). In addition, changes 
in climate can affect the chemistry of ozone, including in- 
creasing the rate of OH production due to higher water va- 
por concentrations with warming temperatures, that will have 
impacts for methane and other gases, e.g., (Stevenson et al., 
2006). Uncertainties in these processes impact knowledge of 
preindustrial ozone concentrations, the evolution of ozone, 
and present-day distributions. These factors contribute to the 
broad range of radiative forcing estimates in the IPCC. 

Ozone can also indirectly affect the climate through the 
carbon and hydrological cycles. As it is a phytotoxin, in- 
creases in ozone can reduce global primary productivity and 
therefore the CO 2 uptake by biota (plants, forests), poten- 
tially leading to an indirect forcing of a similar magnitude 
to the direct ozone forcing (Sitch et al., 2007; Collins et al., 
2010). These interactions have yet to be explored across a 
suite of models. 

Studies since the Intergovernmental Panel on Climate 
Change (IPCC) Fourth Assessment Report (AR4) have 
found that changes in short-lived species, including tropo- 
spheric ozone, are expected to have a significant impact 
on global temperatures, especially over the Northern Hemi- 


sphere (Levy et al., 2008; Shindell et al., 2008). In particular, 
by 2050, Short-Lived Climate Forcing (SLCF) agents could 
be responsible for up to 20 % of simulated global mean an- 
nually averaged warming and up to 40 percent of the total 
projected summertime warming in the central United States 
from 2050-2100 (Levy et al., 2008). However, the contribu- 
tion of short-lived species to simulated changes in global- 
mean surface temperature can have a strong dependence on 
both the SLCF magnitude and spatial distribution (Shindell 
and Faluvegi, 2009). These factors have been the impetus for 
policy formulation that mitigates both air quality and global 
warming, (West et al., 2006, 2007; van Vuuren et al., 2006; 
Wallack and Ramanathan, 2009; Ramanathan and Xu, 2010; 
Shindell et al., 2013) with the concomitant need for accurate 
measurements and model results. 

The Atmospheric Chemistry and Climate Model Inter- 
comparison Project was initiated to complement the Climate 
Model Intercomparison Project, Phase 5 (CMIP5) (Taylor 
et al., 2011) by focusing on the coupling of chemically ac- 
tive atmospheric constituents with climate both historically 
and in the future (Lamarque et al., 2013). A suite of state- 
of-the-art chemistry-climate models were driven by common 
emissions and comparable boundary conditions over differ- 
ent time periods from preindustrial to present-day. In ad- 
dition, simulations of future climate were calculated using 
emissions derived from Representative Concentration Path- 
way (RCP) scenarios (van Vuuren et al., 2011). This ap- 
proach can help quantify the role of SLCF in past and future 
climate response across a range of models, which was diffi- 
cult to do previously (Shindell and Lamarque, 2009). Conse- 
quently, ACCMIP will be a valuable resource for the IPCC 
AR5. 

Furthermore, a constellation of satellites that can observe 
a number of trace gas constituents are now available with 
global coverage and data records of sufficient length to be 
useful for the evaluation of ACCMIP simulations. Expand- 
ing from the analysis of Aghedo et al. (2011b), we evaluate 
how ACCMIP tropospheric ozone bias effects OLR bias us- 
ing global, simultaneous observations of tropospheric ozone 
and sensitivity of OLR to tropospheric ozone from the Tropo- 
spheric Emission Spectrometer, which was launched aboard 
the NASA Aura spacecraft in 2004 (Beer, 2006), using data 
collected from 2005-2010. TES does not measure the full 
OLR but does measure spectral radiances in the 9.6 mi- 
cron ozone band. However, tropospheric ozone has by far the 
strongest infrared absorption in the 9.6 micron band (Clough 
and Iacono, 1995). Consequently, the impact of tropospheric 
ozone on OLR bias can be computed from the 9.6 micron 
ozone band alone. We apply TES Instantaneous Radiative 
Kernels (IRK) to the ACCMIP tropospheric ozone bias in 
order to relate changes in the vertical structure of tropo- 
spheric ozone to changes in OLR (8 OLR) at the top-of-the- 
atmosphere (TOA). The TES IRK quantifies the sensitiv- 
ity of OLR to changes in the vertical distribution of tropo- 
spheric ozone for each satellite observation, accounting for 
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the effects of temperature, water vapor and effective cloud 
properties that are simultaneously estimated (Worden et ah, 
2011). These corresponding state parameters largely deter- 
mine atmospheric opacity and therefore control the magni- 
tude of the IRK. Consequently, we quantify the change in 
OLR due to ACCMIP tropospheric ozone bias referenced 
against the observed atmospheric state as opposed to a model 
simulated atmospheric state. 

We then explore the implications of ACCMIP OLR bi- 
ases from tropospheric ozone on ozone radiative forcing 
(RF), which is defined as the net change in irradiance at 
the tropopause after stratospheric temperature adjustment 
due to changes in ozone concentrations from preindustrial 
to present day (Forster et al., 2007). Ozone RF is depen- 
dent on a number of factors, including preindustrial emis- 
sions, background concentrations (through lightning-derived 
NOx and stratospheric input), present-day concentrations, at- 
mospheric radiative transfer, chemical mechanisms, and the 
physical climate (e.g. clouds, temperature, humidity). Uncer- 
tainties in these processes potentially could lead to biases in 
present day OLR. Hence, we investigate whether ACCMIP 
ozone RF is correlated with ACCMIP 8 OLR. 


2 Tropospheric emission spectrometer 

Launched in July 2004, the NASA EOS Aura platform is in 
a polar, Sun-synchronous orbit with an equator crossing time 
of 13:40 and 02:29 local mean solar time for ascending and 
descending orbit paths, respectively. The Tropospheric Emis- 
sion Spectrometer (TES) is a Fourier Transform Spectrome- 
ter (FTS) that measures spectrally-resolved outgoing long- 
wave radiation of the Earth’s surface and atmosphere at fre- 
quencies between 650-2250 cm -1 with unapodized spectral 
resolution of 0.06 cm -1 (Beer, 2006). This spectral resolu- 
tion is sufficient to resolve pressure-broadened IR absorption 
lines in the troposphere and allows TES to estimate simulta- 
neous vertical profiles of ozone, water vapor, carbon monox- 
ide, methane, deuterated water vapor as well as atmospheric 
temperature. Previous studies used these profiles under clear- 
sky scenes to estimate the greenhouse gas effect from upper 
tropospheric ozone to be 0.48 =b 0.14 Wm -2 (Worden et al., 
2008). Effective cloud pressure and optical depth, surface 
temperature and land emissivity are also derived from TES 
spectral radiances, which allows TES to perform ozone re- 
trievals in all-sky conditions (Eldering et al., 2008). Algo- 
rithms for radiometric calibration (Worden et al., 2007), re- 
trieval of atmospheric parameters (Bowman et al., 2006), 
and error characterization (Worden et al., 2004) along with 
cloud property retrievals (Kulawik et al., 2006) have been 
described previously. The TES forward model used for com- 
puting spectral radiances and Jacobians (Clough et al., 2006) 
is based on LBLRTM (line-by-line radiative transfer model), 
which has been used as the basis for a number of radiative 
transfer models in climate models (Forster et al., 201 1). TES 


radiances have been compared to other satellite and aircraft 
data (Shephard et al., 2008). There is no detectable trend in 
residual spectral radiances (observed minus calculated) from 
2005-2009 to within 0.6 K (Connor et al., 2011). Previous 
validation studies of the TES ozone V002 product indicate 
ozone profiles are biased high in the troposphere (~15%) 
relative to ozonesondes (Nassar et al., 2008) and aircraft data 
(Richards et al., 2008) while total ozone columns are biased 
high by about 10 DU relative to OMI measurements (Oster- 
man et al., 2008). The optimal estimation technique used op- 
erationally on TES data provides formal uncertainties (Bow- 
man et al., 2002, 2006). These uncertainty estimates were 
tested by comparisons with northern latitude ozonesondes re- 
ported in Boxe et al. (2010), which showed that formal mea- 
surement and species interference errors were consistent with 
empirical calculations. 

We use TES V004 ozone profile product for this 
study, which has been updated relative to the V002 prod- 
uct. Ozonesonde profile comparisons with the Interconti- 
nental Chemical Transport Experiment Ozonesonde Net- 
work Study (IONS), the World Ozone and Ultraviolet 
Data Center (WOUDC), the Global Monitoring Division 
of the Earth System Research Laboratory (GMD-ESRL), 
and the Southern Hemisphere Additional Ozonesonde 
archives (SHADOZ) have been performed. Approximately 
5000 matches with one or more of these data sources 
are found using the coincidence criteria of =L3h and a 
300 km radius, spanning a latitude range and time-span 
from 72.5° S to 80.3° N from 2004 to 2008. Analysis 
of the comparisons are documented in the TES Valida- 
tion Report V4 http://eosweb.larc.nasa.gov/PRODOCS/tes/ 
validation/TES _ValidationJReport_v40.pdf. In the Northern 
Hemisphere, the TES tropospheric ozone is biased high by 
about 10-15 % throughout the free troposphere. In the trop- 
ics, the bias in the free troposphere is vertically dependent 
with about a 10 % high bias at about 800 hPa and decreasing 
to near zero at 200 hPa. The southern midlatitudes show sim- 
ilar features but are slightly low biased (<5 %) at 200 hPa. 
Comparisons of TES tropospheric ozone with ozonesondes 
by Young et al. (2013) show a similar sign and magnitude of 
bias. 


3 Instantaneous radiative kernels 


Following Worden et al. (2011), the logarithmic instanta- 
neous radiative kernel (LIRK) is defined as the sensitivity 
of OLR at TOA to changes in the vertical distribution of tro- 
pospheric ozone: 


9^toa 

dlnq(zi) 


2jt 2 

err d Ljq\(v , o , cb) 

J J J ^ — — — - cos OsinO d0 depdv (1) 

V 0 0 


where v is frequency, FYoa is the ozone band radiative flux 
from integrating over v £[985, 1080] cm -1 , q is an ozone 
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Fig. 1 . Logarithmic instantaneous radiative kernels (LIRK) in mWm 2 /ln(vmro 3 ), TES zonal ozone in parts-per-billion (ppb), and longwave 
radiative effect (LWRE) in Wm -2 are shown along across each column respectively. Each row represents data averaged for 2005-2009 for 
December- January-February (DJF), March-April-May (MAM), June- July-August (JJA), and September-October-November (SON). 


profile in volume mixing ratio as a function of altitude z at 
level l , Ltoa is the spectral radiance at frequency v, zenith 
angle 0 and azimuth angle 0. The partial derivatives of spec- 
tral radiance Ltoa are provided by the TES operational ra- 
diative transfer algorithm, (Clough et al., 2006), and are used 
within the retrieval algorithm to estimate the vertical distri- 
butions of trace gases, temperature, water vapor, and clouds 
(Worden et al., 2004; Bowman et al., 2006). Changes in other 
atmospheric parameters are incorporated into the LIRK, e.g., 
retrieved cloud top height will change the altitudes for which 
Ltoa is sensitive to ozone variations. From Eq. (1), increases 
in tropospheric ozone lead to a reduction in OLR consistent 
with the ozone greenhouse gas effect. 

Equation (1) is treated as an operator in a first order Taylor 
series expansion to calculate the ACCMIP OLR bias from 
tropospheric ozone as 

SOLR'i'"' = Hj ( zi ) [In qf ( z, ) - In qf* &)] (2) 

where qj(zi) and q° bs (zi ) are the model and TES ozone, re- 
spectively, at the j th location and the I th altitude level while 


Hj(zi) = 


a pi 

or TOA 

3 In qj(zi) 


( 3 ) 


is the LIRK at location j and evaluated at altitude zi dis- 
cretized at levels /. Consistent with the ozone greenhouse 
gas effect, the LIRK is negative for increases in tropospheric 
ozone so that for a positive tropospheric ozone bias between 
a model and TES at an altitude zu SOLRj ,m < 0. In the sub- 
sequent sections, the level l will also refer to corresponding 
pressure levels. 

The total longwave radiative effect (LWRE) is calculated 
by integrating Eq. (3) in altitude to the tropopause. The 
LWRE can be thought of as the reduction in OLR to a 100 % 
change in the tropospheric ozone profile (Worden et al., 
2011). Note that the LWRE is a logarithmic change refer- 
enced to TES ozone and consequently it can not be used to 
calculate the change to a complete absence of ozone. 

TES ozone, LIRK, and LWRE averaged for 2005-2009 are 
shown in Fig. 1. The LWRE was integrated to the thermal 
tropopause height derived from the Goddard Modeling and 
Assimilation Office (GMAO) GEOS-5 (Molod et al., 2012). 
There is a strong sensitivity in the LIRK of up to 35 m Wm -2 


Atmos. Chem. Phys., 13, 4057-4072, 2013 


www.atmos-chem-phys.net/13/4057/2013/ 




K. W. Bowman et al.: Ozone RF 


4061 


in the mid-troposphere extending as low as 600 hPa and 
arcing poleward to 200 hPa at 60° N and 60° S. The LIRK 
increases in the upper troposphere and lower stratosphere 
with sensitivities well above 60 m Wm" 2 per unity change 
in In q . The LIRK spatial pattern is remarkably similar to cli- 
matological relative humidity (RH) distributions, which are 
at a minimum in subsidence regions especially in the South- 
ern Hemispheric sub-tropics. Interestingly, low RH in the SH 
sub-tropics have been linked to equilibrium climate sensitiv- 
ity (Fasullo and Trenberth, 2012). This similarity suggests 
that the LIRK sensitivity is driven primarily by atmospheric 
opacity which in turn is controlled by a combination of tem- 
perature, water vapor, and clouds. There is a seasonal migra- 
tion of the LIRK maximum across the equator, which follows 
the change in the inter-tropical convergence zone (ITCZ) and 
is driven primarily by the seasonal shift in cloud distribu- 
tions. 

The seasonal pattern of ozone in the SH is strongly in- 
fluenced by the presence of biomass burning and lightning 
leading to a maximum of 60-70 ppb in September-October- 
November (SON). The impact of biomass burning and light- 
ning in South America, sub-equatorial Africa, and Indone- 
sia are clearly seen. Near source regions, the LWRE exceeds 
1 Wm -2 during SON. During December- January-February 
(DJF), ozone from biomass burning in Africa north of the 
(ITCZ) has similarly high LWRE. In most months, there is a 
persistently high LWRE (> 1 Wm -2 ) over the Middle East. 
This is driven in part by relatively clear skies and a strong 
thermal contrast that amplifies the ozone greenhouse gas ef- 
fect. However, in the summer months, there is an ozone en- 
hancement in the middle troposphere (400-500 hPa) induced 
by trapping from Saharan and Arabian anticyclones, which 
also corresponds to the highest magnitude LIRK values in 
the middle troposphere (Li et al., 2001; Liu et al., 2009). 


4 Methodology 
4.1 Radiative forcing 

Radiative forcing (RF) is a measure of the energy imbalance 
of the Earth-atmosphere system and is used as a means of 
quantifying the potential of external agents to perturb that 
system. In the context of historic climate change, the pertur- 
bation is referenced from preindustrial (1750s) to the present- 
day (2000s). In order for the ozone forcing to be taken as 
a reasonable proxy for the expected temperature response, 
RF is specifically defined as the change in net irradiance 
at the tropopause after stratospheric temperatures have re- 
laxed to radiative-dynamical equilibrium but with the sur- 
face and atmospheric state held fixed (Forster et al., 1997, 
2007). The forcing-response relationship, however, can be 
more complicated depending on the forcing agent and its 
distribution (Fels et al., 1980; Hansen et al., 1997, 2007; 
Shindell and Faluvegi, 2009). Methods of calculating tropo- 


spheric ozone RF from models using “off-line” techniques, 
i.e., radiative transfer calculations performed independently 
of a climate model’s internal radiation calculation, follow 
two steps. The first step is to calculate the change in ozone 
concentrations by forcing a global chemistry-climate model 
with pre-industrial (usually taken in the 1750s) and present- 
day concentrations in separate “time-slices”, which is an in- 
terval of time, e.g., 1750-1760. The second step is to calcu- 
late the change in longwave and shortwave radiation due to 
the change in present-day ozone relative to the pre-industrial 
era using an “off-line” radiative transfer model (RTM) where 
stratospheric temperatures are allowed to equilibrate result- 
ing in an adjusted irradiance (Edwards and Slingo, 1996; 
Stevenson et al., 2006; Knutti and Hegerl, 2008; Steven- 
son et al., 2013). Radiative forcing is generally referenced 
at the tropopause, the definition of which can have a sig- 
nificant impact on the final calculation, e.g. flat tropopause 
set at 100, 150 or 200 hPa, zonally invariant and linear with 
latitude tropopause (Naik et al., 2005; Hansen et al., 2007), 
chemical tropopause using the 150 ppbv ozone level (Steven- 
son et al., 2006), and the WMO thermal tropopause (Aghedo 
et al., 201 lb). Radiative forcing is simply the difference be- 
tween the irradiance at two different time periods, which can 
be expressed as: 

RF m = R ( qp - R(0 (4) 

where R is the globally, area-weighted average net irradi- 
ance (shortwave (S W) plus longwave (LW)) in m Wm -2 in- 
cluding stratospheric readjustment, is the vertical distri- 
bution (discretized as a vector) of present-day ozone simu- 
lated by model m, and q™ is preindustrial ozone for the same 
model. Implicit in Eq. (4) is temporal averaging over some 
time-slice. 

4.2 Application to TES data 

TES observations directly measure the 9.6 micron ozone 
band radiances (the primary band where ozone absorbs ther- 
mal infrared radiation) and have the spectral resolution to 
disentangle the geophysical quantities, e.g., temperature and 
clouds, driving ozone band radiative flux variability. These 
data have the potential to assess model results with respect to 
the longwave component of ozone radiative forcing. How- 
ever, differences between satellite and climate model cal- 
culations need to be considered. TES only measures ther- 
mal radiance at the top-of-the-atmosphere (TOA). Conse- 
quently, these measurements do not provide information on 
the shortwave irradiance. Chemistry-climate models calcu- 
late the full global and diurnal cycle of ozone whereas TES 
has a global repeat cycle of 16 days and can only measure 
twice a day through its ascending and descending nodes. 
However, Aghedo et al. (2011a) showed that TES sampling 
is sufficient to capture zonal scale variations in ozone when 
compared to chemistry-climate model simulations. 
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The mean ACCMIP bias in OLR from tropospheric ozone 
is calculated from Eq. (2) as 

SOLRi = — Y Y2 w i - ln[q° bs ],- > /) (5) 

ieDj leL 

where Wi are area-weights (to account for the relative areas 
of different latitude bands), Dj is a set of observed locations, 
Nj is the number of locations in set Dj , and L is the set 
of altitude levels whose maximum value is at the tropopause, 
which we choose to be the chemical tropopause q = 150 ppb. 
The domain Dj can be comprised of either grid cells or zonal 
bands. In either case, we will define the OLR bias as <50LR4 
where j refers to latitudinal bands or individual grid boxes 
depending on the context described in the subsequent sec- 
tions. The global mean OLR bias from tropospheric ozone is 
calculated by averaging over all Dj and will be denoted as 
<50LR m . The ACCMIP ensemble mean OLR bias from tro- 
pospheric ozone will be denoted as 5 OLR j for locations Dj 
and the ensemble global mean is simply 8 OLR. Implicit in 
Eq. (5) is temporal averaging from 2005-2010. 


5 Results 

5.1 ACCMIP simulations 

We use TES observations to evaluate the OLR bias from 
tropospheric ozone in the chemistry-climate models that 
participated in ACCMIP. A complete description of the 
chemistry-climate and chemical transport models along 
with their short hand designation (CESM-CAM-superfast, 
CICERO-OsloCTM2, CMAM, GEOSCCM, GFDL-AM3, 
GISS-E2-R, GISS-E2-TOMAS, HadGEM2, LMDzOR- 
INCA, MIROC-CHEM, MOCAGE, NCAR-CAM3.5, 
STOC-HadAM3, UM-CAM) can be found in Lamarque 
et al. (2013). Each were driven by a common set of emis- 
sions (Lamarque et al., 2010). The experimental design 
was based on decadal “time-slice” experiments driven by 
decadal mean sea surface temperatures (SST). The historic 
periods included 1850, 1980, and 2000 along with 2030 and 
2100 time slices that follow the representative concentration 
pathways (RCP) (van Vuuren et al., 2011). As described 
by Lamarque et al. (2013) the level of complexity in the 
chemistry schemes varied significantly between the models. 
The physical climate was based on prescribed SSTs for 
most models with the notable exception of the GISS model, 
which is integrated with a fully coupled ocean-atmosphere 
model (Shindell et al., 2013). Following Young et al. (2013), 
model simulations for the 2000 decade were averaged and 
then interpolated to the domain of the gridded TES product 
archived along with the ACCMIP simulations at the British 
Atmospheric Data Archive (BADC), which is composed 
of 64 pressure levels at 2 x 2.5° spatial resolution. These 
decadal mean model simulations are then compared against 
TES tropospheric ozone for 2005-2010. 
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Fig. 2. Zonal bias between ACCMIP models and TES ozone aver- 
aged over 2005-2010. Ensemble average of ACCMIP compared to 
TES ozone is shown under ENS in the bottom right. 


5.2 Application to ACCMIP models 

We first compare the zonal-vertical difference between the 
ACCMIP and TES tropospheric ozone from 2005-2010 as 
shown in Fig. 2 including the ensemble (ENS) in the bottom 
right. There is a rich diversity in the zonal ozone distribution. 
In the middle to lower troposphere the agreement is gener- 
ally within 10-15 ppb. In the Northern Midlatitudes (NMLT), 
the GISS-E2-R, GISS-E2-TOMAS, and MOCAGE were 
larger than TES ozone estimates by more than 10 ppb. On 
the other hand, CICERO-OsloCTM2, HadGEM2, MIROC- 
CHEM, and CMAM tend to underestimate NMLT ozone rel- 
ative to TES ozone. In the tropical troposphere, most mod- 
els tend to underestimate ozone. The upper troposphere and 
lower stratosphere show stronger differences. The CMAM 
and MIROC-CHEM models have significantly higher ozone 
in both the NMLT and the Southern Midlatitudes (SMLT). 
Conversely, MOCAGE, HadGEM2, STOC-HadAM3, and 
UM-CAM estimate lower ozone than TES ozone in both the 
NMLT and SMLT. However, the TES observation operator is 
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Fig. 3. Zonal distribution of ACCMIP OLR bias from tropospheric ozone, <$OLR^’ m , from 2005-2010. The vertical scale is defined from 
1000-100 hPa discretized to TES pressure levels. Color scale is inverted to show that positive <50LR^ ,m is a model low ozone bias relative to 
TES. Ensemble mean bias <$0LR/ is denoted by ENS. 


not applied to these models as Aghedo et al. (201 lb), nor is 
there any bias correction applied to the TES data. Based on 
ozonesonde comparisons discussed in Sect. 2, TES ozone is 
systematically high in the Northern Hemisphere, which sug- 
gests that the high ozone biases in GISS-E2-R and MOCAGE 
are exacerbated but the slight low ozone bias in the ensemble 
distribution is likely insignificant. TES upper tropospheric 
ozone in the tropics is consistent with ozonesonde measure- 
ments, which indicates that the ensemble model low bias is 
robust. 

The ACCMIP OLR bias from tropospheric ozone, 
<$OLR/’ m , which is calculated from Eq. (2), is shown in Fig. 3 
discretized at TES pressure levels. The zonal cross-section in 
this figure can be related to the ACCMIP ozone bias in Fig. 2 


through Eq. (2). For example, in Fig. 2 a lOppb low bias 
of GISS-E2-R relative to TES ozone at 300 hPa and 15° S 
leads to a 80LRj ,m of about 5mWm -2 at the same loca- 
tion as shown in Fig. 3. As quantified by the LIRK in Fig. 1, 
the thermal contrast between the temperature at which the 
ozone absorbs in the troposphere and the surface temper- 
ature along with the low humidity, which decreases atmo- 
spheric opacity, contributes to the importance of the tropics 
and subtropics relative to the high latitudes. Consequently, 
underestimates of tropical and subtropical ozone are ampli- 
fied in terms of 80LRj ,m relative to the extratropics. Almost 
half of the ACCMIP models have 80LRj ,m that exceeds 
lOmWm -2 at individual pressure levels /. Radiatively sig- 
nificant differences are not confined to the upper troposphere. 
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Fig. 4. Spatial distribution of ACCMIP OLR bias from tropospheric ozone, SOLR^, from 2005-2010, limited to 80° S-80° N, and based on 
a chemical tropopause q = 150ppb as diagnosed from TES. ACCMIP ensemble, SOLRj , is denoted by ENS. 


In several models, tropical differences in ozone at pressures 
greater than 600 hPa lead to 80LRj ,m > lOmWm -2 . 

The ACCMIP OLR bias from tropospheric ozone is the 
ACCMIP tropospheric ozone bias weighted by the LIRK. 
Comparing the LIRK distributions in Fig. 1 to the 80LRj ,m 
in Fig. 3 suggest that in the tropics and SH, <50LR^’ m from 
models such CESM-superfast, HadGEM-2, MOCAGE, UM- 
CAM, NCAR-CAM 3.5, STOC-HadAM3 has a strong in- 
fluence from the LIRK, particularly the positive 50LR j' m 
at 500 hPa around 10-15°S and extending poleward up to 
100 hPa around 50° S. This LIRK spatial pattern is persistent 
across a number of models in SH as can be seen in the AC- 
CMIP ensemble mean <50LR/ in Fig. 3. In the NH high lati- 
tudes, such as with GISS-E2-R at 60° N and 600 hPa, there is 
a significant bias in ozone but the attenuated OLR sensitivity 
makes SOLR.j’ m relatively small. 

A complimentary perspective is shown in Fig. 4, which 
shows <$OLr 4 calculated from Eq. (5) where <$OLR/’ m at 
each location j is vertically integrated up to the chemical 


tropopause of q = 150 ppb diagnosed from TES ozone. Most 
ACCMIP models have a positive 50LR4 in the SH trop- 
ics because they underestimate tropospheric ozone relative 
to TES tropospheric ozone. The overestimate in the Eastern 
Tropical Atlantic is a persistent feature in all of the ACCMIP 
models and is reflected in the ensemble mean. Tropospheric 
ozone distributions in the tropical Atlantic are driven by a 
number of processes but are dominated by lightning (Jacob 
et al., 1996; Jenkins and Ryu, 2004; Sauvage et al., 2007; 
Bowman et al., 2009). A second persistent feature is the 
overestimate centered over Southern Africa. Southern equa- 
torial Africa is an important source of biomass burning (Ed- 
wards et al., 2006; Aghedo et al., 2007). Satellite-based “top- 
down” estimates indicate the emissions from biomass burn- 
ing are significantly underestimated (Arellano et al., 2006; 
Jones et al., 2009). Local sources of pollution and biomass 
burning have been associated with upward trends in ozone, 
particularly in the lower troposphere (Clain et al., 2009). 
However, the Southern African region has a complex circu- 
lation pattern that includes both anticyclonic transport and 
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Fig. 5. Zonal SOLR^ between ACCMIP and TES from 2005-2010, 
limited to 80° S-80° N, and is based on a TES diagnosed chem- 
ical tropopause ( q = 150ppb). The “ENS” refers to the ACCMIP 
ensemble average, <$OLRy. Positive values indicate model ozone is 
biased low relative to TES ozone and consequently the model OLR 
is biased high relative to TES OLR. 


recirculation as well as direct eastward and westward trans- 
port (Garstang et al., 1996; Sinha et al., 2004). Upwind 
sources of ozone precursors from biomass burning, pollu- 
tion, and lightning can be advected across Southern Africa 
and out to the remote Pacific (Chatfield and Delany, 1990; 
Chatfield et al., 2002). Some of the radiatively strongest dif- 
ferences (>200 mW/m 2 ) in <50LR4, e.g, GISS-E2-R and 
NCAR-CAM3.5, are throughout the tropical Pacific, which 
may contribute to the positive SH <50LR ; . The ACCMIP 
ensemble shows a persistent pattern in the tropical Pacific 
with a local maximum near 75 m Wm -2 . Comparison of the 
tropical distribution in Fig. 4 with the zonal-vertical distribu- 
tion in Fig. 3 of SOLR^, points to the combination of low 
ozone throughout the troposphere that is amplified by the 
strong mid-tropospheric radiative sensitivity of the southern 
branch of the LIRK in Fig. 1. Tropical ozone is sensitive to 
convective mass flux, height, and subsidence, particularly in 
the Eastern Pacific (Liu et al., 2010). These factors in con- 
junction with ozone precursor uncertainties could help ex- 
plain these features, though they are more prevalent during 
El Nino periods such as in 2006 (Nassar et al., 2009; Chan- 
dra et al., 2009) which are not well simulated in ACCMIP 
due to decadally averaged SST boundary conditions. 

The SH tropics and subtropics dominate radiative dif- 
ferences as shown by the vertically integrated zonal dis- 
tribution in Fig. 5. The ensemble mean <50LR ; is high 
by about 120 mWm -2 between 10-20° S. The mean 
bias is strongly influenced by MOCAGE, CESM-CAM, 
and NCAR-CAM3.5, whose (SOLR^ varies from 200- 
400 mWm -2 though the latter two are driven by similar 
physical climate models. In the NH extratropics, the en- 
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Fig. 6. ACCMIP ozone radiative forcing from preindustrial to 
present day as a function of SOLR m , which is the difference be- 
tween ACCMIP model ozone and TES ozone weighted by the TES 
Instantaneous Radiative Kernel. The 50LR is integrated through 
the TES diagnosed chemical tropopause q = 150ppb and averaged 
from 2005-2010. 


semble mean SOLRj ~ OmWm" 2 with the extrema near 
lOOmWm -2 at 30° N decreasing poleward at a rate be- 
tween linear and exponential. Systematic overestimates in 
the northern high latitudes in the TES ozone data will not 
strongly effect the ensemble mean bias because of the lower 
OLR sensitivity to tropospheric ozone. On the other hand, the 
SH low bias is robust because of the small upper tropospheric 
system errors in the TES retrieval. 

5.3 Implications for ACCMIP ozone radiative forcing 

We explore whether present-day OLR model bias from tro- 
pospheric ozone is correlated with model ozone RF and 
whether this relationship can be used to reduce our uncer- 
tainty in preindustrial to present day ozone RF. The cause of 
ACCMIP OLR bias from tropospheric ozone is a complex 
combination of uncertainties related to both the physical and 
chemical climate. Some of these differences are related to 
natural background conditions, e.g., lightning NO x . Model 
biases for present-day background ozone may be expected to 
be correlated with respective biases from pre-industrial back- 
ground ozone. Young et al. (2013) showed, for example, that 
ACCMIP models with high ozone burdens over the present- 
day had high burdens for the other periods of time including 
the preindustrial period. Uncertainties in present day emis- 
sions and their chemical transformation could result in sig- 
nificant biases in tropospheric ozone. Consequently, present 
day measurements have the potential to inform preindustrial 
to present day ozone radiative forcing. 

In Table 1, the domain Dj of SOLR J m in Eq. (5) is de- 
fined for the tropics (15° S-15° N), NH and SH extra-tropics. 
Tropical biases contributed the most to the global <SOLR m 
with some zonal (SOLR^ >250mWm -2 . There are 5 of 
the 15 models with tropical <50LR^ > lOOmWm -2 . While 
most of the models underestimate ozone in the tropics-and 
hence have a positive 5 OLR J m , there are 2 models that have 
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Table 1. ACCMIP OLR bias from tropospheric ozone (<$OLR;4) and ozone radiative forcing (RF m ) in m Wm -2 . The ozone RF m is taken 
from (Stevenson et al., 2013) and is defined relative to 1750 and at a climatological thermal tropopause. The ACCMIP ensemble mean and 
1-cr are shown in the last column of the first two rows. Tropics are defined from 15° S-15° N. 


Model 

Tropics <$OFr4 

SH SOLR J m 

NH SOLR'Ih 

Total 50FR m 

RF m 

ENS 

79 

47 

12 

39 

389 

ENS STDev 

75 

49 

32 

41 

60 

CESM-CAM-superfast 

256 

103 

59 

115 

468 

CICERO-OsloCTM2 

81 

49 

42 

52 

423 

CMAM 

106 

11 

23 

35 

355 

GEOSCCM 

29 

2 

-16 

0 

404 

GFDF-AM3 

1 

-8 

-18 

-10 

446 

GISS-E2-R 

67 

24 

-41 

6 

326 

GISS-E2-R-TOMAS 

65 

6 

-37 

0 

345 

HadGEM2 

50 

75 

33 

53 

341 

FMDzORINCA 

-6 

13 

-6 

1 

384 

MIROC-CHEM 

-37 

20 

16 

7 

416 

MOCAGE 

129 

169 

48 

112 

250 

NCAR-CAM3.5 

168 

54 

26 

65 

446 

STOC-HadAM3 

101 

85 

37 

69 

436 

UM-CAM 

96 

61 

-1 

43 

411 



Fig. 7. The deviation of ACCMIP ozone radiative forcing from an 
ensemble mean versus |50LR m |, which is the magnitude of the 
difference between model ozone and TES ozone weighted by the 
TES Instantaneous Radiative Kernel. The mean is calculate from 
an subset of models whose |<$OLR m | < 10m Wm -2 . The <$OLR 
is integrated through the TES diagnosed chemical tropopause q = 
150ppb and averaged from 2005-2010. The correspondence be- 
tween the colors and models is the same as in Fig. 6. The brackets 
for the slope represent the uncertainty at the 95 % confidence. 


a negative (SOLR^. Similarly, the SH SOLRh reaches as 
much as 170mWm~ 2 . The NH has a much wider response 
with about half of the models having a SOLR^ as high as 
60 m Wm -2 . There are 5 models whose total <50LR m are less 
than 10 m Wm -2 whereas most are greater than 30 m Wm -2 . 
The ensemble mean and standard deviation for <50LR m is 
8 OLR = 39 =b 41 m Wm -2 , which reflects this variability. 

In order to explore the implications of these biases 
on ozone radiative forcing, the global <50LR m is com- 
pared to ACCMIP ozone RF in Fig. 6. The ACCMIP 


ozone RF is based on the climatological thermal tropopause 
and the ozone RF has been extended from 1850 to 
1750 by adding 40mWm -2 (Skeie et al., 2011; Steven- 
son et al., 2013). There is a wide range of ozone 
RF estimates from 210-428 m Wm -2 with a mean of 
384 =b 60m Wm -2 . To investigate whether there is a re- 
lationship between model ozone RF and the magnitude 
of 50LR m we first choose a subset of models whose 
|<$OLR m | < 10m Wm" 2 (GEOSCCM,GISS-E2R,GISS-E2- 
R-TOMAS,LMDzORINCA,MIROC-CHEM). The ozone 
RF mean and standard deviation of this subset is 
375 zb38.4mWm -2 . We then calculate the ozone RF abso- 
lute deviation of all the ACCMIP models with respect to the 
subset mean as 

A m = \R s -R m \ (6) 

where R s is the subset ACCMIP ensemble mean ozone RF. 
The ACCMIP deviation is shown as a function of |50LR m | 
in Fig. 7. The color scheme in Fig. 7 is the same as in 
Fig. 6. The slope of the linear fit shown in Fig. 7 is 0.60 
with a 95% confidence interval of [0.29,0.91] and R 2 = 
0.59. This correlation is largely driven by two models whose 
|<50LR m | > lOOmWm -2 . If these models are not included 
then R 2 = 0.18 and the null hypothesis, i.e., there is no sig- 
nificant relationship between ACCMIP OLR bias and ozone 
RF, can only be rejected at the 80% confidence level. It’s 
interesting to note that these two models - CESM-CAM and 
MOCAGE - have the highest (538) and lowest (309), respec- 
tively, ozone radiative forcing. The ozone RF is driven in part 
by a change of only 4.8 Dobson Units (DU) in ozone for 
MOCAGE but a 10 DU ozone change for CESM-CAM from 
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preindustrial to present-day (Stevenson et al., 2013; Young 
et al., 2013). 

The difference in correlation between the ACCMIP en- 
semble subgroups provides a reasonable basis for ensemble 
selection. As a limitation for ensemble selection, TES obser- 
vations can not be used to distinguish between models for 
which there is no correlation between model OLR bias from 
tropospheric ozone and model ozone RF. Therefore, models 
are selected for which |50LR m | < lOOmWm -2 and there- 
fore A m < 86mWm -2 . The ACCMIP radiative forcing un- 
der these conditions is 394 ±42 m Wm" 2 . The mean is about 
the same and the standard deviation has been reduced by 
28% relative to 389± 60mWm -2 , which is the ACCMIP 
ensemble mean ozone RF in Table 1. While this estimate 
has reduced variability relative to the full ACCMIP ensem- 
ble ozone radiative forcing, it does not fully account for other 
sources of uncertainty, e.g., differences in radiative transfer 
models (Forster et al., 2011). 

6 Conclusions 

We have presented an evaluation of ACCMIP OLR from tro- 
pospheric ozone using TES ozone and LIRK data. We fur- 
ther explored the implication of this evaluation on ozone 
RF estimates. We find a significant correlation R 2 = 0.59 
between |<$OLR m | and the absolute deviation from the AC- 
CMIP ensemble ozone RF, A m . However, this correlation 
drops to R 2 = 0.18 if we only include models for which 
|<50LR| < 100 mWm -2 . Using this change in correlation as 
a basis for ensemble selection, we estimate the ACCMIP RF 
to be 394±42 m Wm -2 (1 standard deviation), which is close 
to the ACCMIP full ensemble mean RF of Stevenson et al. 
(2013) but with about 30 % less standard deviation. We note 
that our calculation used only 14 of the 16 models used by 
Stevenson et al. (2013). This result is in contrast to the lack 
of correlation across the full ACCMIP ensemble between 
present day tropospheric ozone bias and the change in ozone 
burden from preindustrial to present day reported by Young 
et al. (2013). However, our estimate does not account for 
bias associated with preindustrial emissions nor with others 
sources of uncertainty such as differences in the physical cli- 
mate or radiative transfer algorithms. The results further sug- 
gest that any causal relationship between present day model 
OLR bias and model ozone RF could be model specific and 
not necessarily reflective of any general relationship between 
present day OLR and ozone RF. 

We investigated the spatial patterns driving present-day 
bias in ACCMIP OLR from tropospheric ozone. These pat- 
terns were driven in the ACCMIP models by variations in 
the SH tropics and sub-tropics where the ensemble mean 
5 OLR j reached lOOmWm -2 in some regions. Persistent 
patterns of 50LR^ were centered over the tropical Atlantic 
and over Southern Africa. While the importance of upper tro- 
pospheric ozone to ozone RF is known (Lacis et al., 1990; 


Gauss et al., 2003), significant <50LR j’ m (> lOmWm 2 ) 
could be attributed to ACCMIP ozone bias at pressures ex- 
ceeding 600 hPa. While the modest low bias in 50LR^ in 
the NH can be attributed in part to a systematic TES retrieval 
overestimate, the ACCMIP SH low bias in 5 OLR J m is robust 
in light of the smaller TES tropical retrieval systematic er- 
rors. 

Our results demonstrate a motivation for developing more 
sophisticated approaches that incorporate explicit statistical 
modeling of present day OLR bias from tropospheric ozone 
and ozone RF. Such approaches would permit a more rigor- 
ous means of weighting model ozone RF based upon their 
relationship to observations (Tebaldi et al., 2005; Berliner 
and Kim, 2008). There is considerable interest in using both 
data and an ensemble of climate models to understand his- 
toric change and probabilistically weight future projections 
(Knutti et al., 2002; Collins, 2007; Tebaldi and Knutti, 2007). 
There has been little to no application of these methodologies 
to chemistry-climate projections or the attribution of historic 
climate change to chemically active agents as has been done 
in detection and attribution studies for temperature (Hegerl 
et al., 1996; Santer et al., 2007; Huber and Knutti, 2012). 
There is also the potential of using advanced data assimi- 
lation techniques to attribute observed OLR to natural and 
anthropogenic emissions at finer spatial resolution (Bowman 
and Henze, 2012). The combination of these approaches with 
a process-based analysis (Eyring et al., 2005, 2006; Waugh 
and Eyring, 2008) can help test radiatively important pro- 
cesses against observations in a manner that can reduce our 
uncertainty in ozone RF and increase the reliability of future 
projections. 
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